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Abstract 

In the previous paper, the number distribution of the low-lying spectra around 
Gaussian solutions representing various dimensional fuzzy tori of a tensor model 
was numerically shown to be in accordance with the general relativity on tori. 
In this paper, I perform more detailed numerical analysis of the properties of 
the modes for two-dimensional fuzzy tori, and obtain conclusive evidences for the 
agreement. Under a proposed correspondence between the rank-three tensor in 
tensor models and the metric tensor in the general relativity, conclusive agreement 
is obtained between the profiles of the low-lying modes in a tensor model and 
the metric modes transverse to the general coordinate transformation. Moreover, 
the low-lying modes are shown to be well on a massless trajectory with quartic 
momentum dependence in the tensor model. This is in agreement with that the 
lowest momentum dependence of metric fluctuations in the general relativity will 
come from the i?^-term, since the i?-term is topological in two dimensions. These 
evidences support the idea that the low-lying low-momentum dynamics around the 
Gaussian solutions of tensor models is described by the general relativity. I also 
propose a renormalization procedure for tensor models. A classical application of 
the procedure makes the patterns of the low-lying spectra drastically clearer, and 
suggests also the existence of massive trajectories. 
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1 Introduction 



Thought experiments in combination of quantum mechanics and general relativity [I]- [7] and 
also in string theory [S] show various bounds on accuracy of space-time measurements. The 
existence of such bounds implies that space-time can never be observed as a smooth continuous 
manifold, which is the classical space-time notion in general relativity. Moreover the obvious 
distinction between space-time and matter fields on it would be an obstacle in pursuit of 
unification. Especially, it would be nice if space-time and metric tensor can merely be regarded 
as different sides of a single object. Thus search for a satisfactory alternative to the classical 
space-time notion has theoretical interests. 

The notion of fuzzy space would be a possible candidate for such an alternative [9]- [25]. A 
fuzzy space is defined by the algebra of functions on it, which is allowed to be not only non- 
commutative but also nonassociative. Therefore a fuzzy space is a generalized notion which 
includes both the noncommutative [9]-[l7j and the nonassociative |l8]-[25] spaces. Since an 
algebra can be defined by a rank-three tensoiQ, dynamical models of a rank-three tensor, which 
are denoted by tensor models [26]- [33] in this paper, may be regarded as dynamical theory of 
fuzzy spaces |3l] . An important advantage of the formalism is that fuzzy analogue of the gen- 
eral coordinate transformation can easily be embedded into such tensor models. Another one 
is that one can treat various types of spaces in a unified manner, e.g. irrespective of topologies 
and dimensions. This possibility to describe dynamical fuzzy spaces in terms of tensor models 
has been pursued by the present author, and some successful results have been obtained. In 
[Ml [35l [36l [37] , various classical solutions of tensor models representing fuzzy spaces with some 
physical interests have been obtained, and their properties have been analyzed. Especially in 
the previous paper [27], the number counting of the low- lying fluctuation modes around a 
specific class of solutions, which represent various dimensional fuzzy flat tori and are denoted 
by Gaussian solutions in this paper, agrees with what is expected from the general relativity. 
This result suggests that the effective low-lying physics around the Gaussian solutions can be 
described by the general relativity. If this truly holds, metric tensor and space come from a 
single object as mentioned in the preceding paragraph, because the rank-three tensor is the 
only dynamical variable of tensor models. 

Although the agreement is non-trivially realized in various dimensions, it is obvious that 
such number counting is not enough for a definite conclusion. The main purpose of the present 
paper is to go beyond the number counting to conclusively show the agreement. The paper 
is organized as follows. In Section [21 fuzzy space, tensor models, and Gaussian solutions are 
reviewed. A tensor model with Gaussian solutions, which will be analyzed throughout the 
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present paper, is given. In Section [31 a proposal is given of a correspondence between the 
rank-three tensor of tensor models and the metric tensor in the general relativity. Then the 
DeWitt supermetric with its unique parameter being fixed is derived from a quadratic invariant 
metric of tensor models. The supermetric determines the metric fluctuations transverse to the 
general coordinate transformation. In Section HI numerical analysis of the fluctuation modes 
is performed for two-dimensional fuzzy flat tori. Conclusive agreement is obtained between 
the transverse metric fluctuations and the low-lying fluctuation modes of the tensor model. 
The low-lying spectra are observed to form a massless trajectory with quartic momentum 
dependence. In Section |5l a renormalization procedure for tensor models is proposed. A 
classical application is shown to make much clearer the patterns of the low-lying spectra, and 
also shows the existence of a massive trajectory well over the massless one. The final section 
is devoted to summary, conclusions, and discussions. 

2 Tensor models 

2.1 Fuzzy spaces and tensor models 

The idea of regarding tensor models as dynamical theory of fuzzy spaces was originally pre- 
sented in [34] . In this subsection, I will recapitulate the discussions, stressing more on physical 
motivations. 

Before describing a fuzzy space, let me start with a usual space with coordinates, (fi = 
1, ■ • • ,D). Let me only consider the real functions on the space. One will be able to take a 
basis of the space of all the real functions as {fa{x)} with an appropriate index set. In some 
cases, one can choose a basis with indices having some physical meanings such as frequencies, 
but in general the indices are merely abstract labels. The algebraic relations among the 
functions can be parameterized by a real rank-three tensor Cab'^ as 




(1) 



Introducing a rank-two tensor gab defined by 




(2) 



one obtains 




(3) 



where C, 



abc 
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Here both gab and Cabc are real symmetric tensors, because of the associativity and com- 
mutativity of the products among the functions on a usual space. It is important to note that 
one can freely take another basis for the function space. Since such a change of basis can be 
described by a general linear transformation, tensors related by a general linear transformation 
describe an identical space. The general coordinate transformation can be embedded into the 
general linear transformation, because any function can be expressed as a linear combination 
of the basis functions. Namely, there always exist Mj' such that faix') = Mj'fh{x), where 
x' = x'{x) is a transformed coordinate and Ma^ are the numerical coefficients of expansion. 

If Qab and Cabc are given and gab is invertible, one can get Cab'^ in the product rule ([T]). 
As will be explained later at the end of this subsection, it is more convenient to consider the 
two tensors, g"''' = {g~^)"''^ and Cabc, as the fundamental variables of tensor models rather than 
Cab'^ itself in the product rule. The main reason comes from the requirement of the invariance 
of an action under the general linear transformation, as will be explained at the end of this 
subsection. 

In field theory, the basic variables are fields, which are similar to functions. Their product 
and integration over a space-time are basic operations necessary in extracting physics from 
field theory, so that the tensors g"''^ and Cabc are more directly related to observables in field 
theory than a coordinate system itself. Therefore it would be acceptable in physics to describe 
a space in terms of these tensors, but not by coordinates. 

A fuzzy space may be defined by a deformation of the above tensors away from the values 
corresponding to a usual space. The type of deformation relevant in this paper is basically 
a kind of truncation of the function space by introducing a cutoff to high frequency modes 
without changing the symmetric properties of the tensors g""^ and Cabc under exchange of the 
indices. The cutoff may also be introduced by a cutoff function which vanishes gradually at 
high frequencies. Such truncation of high frequency modes will make it hard to distinguish 
nearby points, which makes a space to become "fuzzy" . 

Because the symmetry under exchange of the indices of the tensors is assumed in the 
deformation, the product rule ([T]) remains commutative. However, this kind of deformation 
generally introduces nonassociativity into the algebra ([1]), because truncation of modes breaks 
closure of an original algebra. The physical importance of this kind of nonassociative fuzzy 
spaces comes from the fact that metric can be incorporated into the function algebra, as will 
explicitly be presented in Section [31 This is in sharp contrast with a noncommutative space 
[9]-[T7j, where metric is given by an additional element, Laplacian [12j . 

The above discussions naturally lead one to suspect that a dynamical theory of fuzzy 
spaces may be constructed as a dynamical theory of the symmetric tensors g"''' and Cabc, and 
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may reproduce the general relativity in a certain limit. Since a fuzzy space is invariant under 
the transformation of a basis, tensor models must have a symmetry under the general linear 
t r ansf ormat ion , 

g^' = {M-')/{M~\^g'^''\ (4) 

which is a fuzzy analogue of the general coordinate transformation on a usual space. Therefore 
an action of a tensor model is a function of g'^^ and Cede-, 

Sig^'^C^de), (5) 

where the upper and lower indices are contracted for the symmetry requirement. In principle, 
it is possible to take another set of dynamical variables than g""^ and Cahc- However the above 
choice seems to be the simplest. For example, one would try to take Cab' as its only dynamical 
variable. Then to construct an invariant, one would need a tensor with more upper indices 
than lower ones. One would be able to make such a tensor by an inverse of a matrix Cah^Ca'c, 
but this would lead to an action singular at the non-invertible values. 

2.2 Euclidean models 

The general discussions on tensor models in the previous subsection seem to favor the existence 
of two kinds of symmetric tensors (^"^ and Cabc as dynamical variables in tensor models. 
However, in the actual numerical analysis, treating both as dynamical variables would become 
too complicated. Also from theoretical viewpoints, the existence of two independent tensors 
would introduce more ambiguities into the formalism. Therefore a Euclidean tensor model 
[35l [37] , which has the non-dynamical rank- two tensor fixed at g""^ = 5°'', will be considered 
throughout this paper. 

It is important to note that such a model is only invariant under an orthogonal subgroup of 
the general linear transformation. Therefore the direct link between the symmetry of a tensor 
model and the general coordinate transformation is missing in such a model. However there 
seem to exist several reasons for the possibility that one may ignore g"'^ as dynamical variables 
without changing the essential features of the system as listed in the following. Especially it is 
noteworthy that, because of the following first reason, Euclidean tensor models are, so called, 
background independent theory of space. 

• Under the assumption that 

gab 

positive definite matrix, it can be diagonalized 
to g"''' = 5""^ by the general linear transformation. This is a partial gauge fixing of the 
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general linear transformation to the orthogonal subgroup. Therefore, classically, there 
are no distinctions between a full tensor model and a Euclidean one under the positivity 
assumptioijl]. In fact, the correspondence to the general relativity fl2T|) . which will be 
discussed later, does not require g"''' to be varied. 

• In the previous paper [37] , it is numerically shown that the distribution of the low-lying 
modes around the Gaussian classical solutions of a Euclidean tensor model agrees with 
the general relativity. 

• What is really required in a tensor model would be a symmetry which resembles the 
general coordinate transformation at low-momentum modes. In general, there are no 
clear distinctions between the general linear and the orthogonal transformations, if the 
two are compared only at a small window of low-momentum modes. 

• The number of components of g"''' is negligible in comparison with Cabc in the limit n — »• 
oo, where n is the total number of the possible values of an index, e.g. a = 1, 2, ■ ■ ■ , n. 
Therefore g""^ would not play essential roles in the dynamics in the limit. 

Thus, in the following discussions, only a Euclidean tensor model, 

S{g^' = 5ab,Cabc), (6) 

will be considered. 

2.3 A Euclidean tensor model with Gaussian classical solutions 

The orthogonal group symmetry of Euclidean models is so loose that the explicit form of 
an invariant action ([6]) has infinite possibilities. In general the dynamics of the models will 
heavily depend on its specific form, and the models have no powers to predict any quantum 
gravitational phenomena. Presently no principles are known to constrain the choices, but 
better understanding of the properties of the models might finally lead to a hint for a preferred 
formulation. It might also be possible that the models can be classified into a finite number 
of universality classes, when the thermodynamics or quantum properties are investigated. 

In this regard, an obvious direction of study is to consider the simplest choices of the 
actions. In the papers [Ml ESI ES], some actions with quadratic and quartic terms in Cabc 
are considered. Especially in [36], considered is an action the equation of motion of which 

^Quantum mechanically, however, one would need to take into account the contributions from the FP ghosts 
associated with the partial gauge fixing. 
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has various commutative but nonassociative deformations of usual spaces, tori and spheres of 
various dimensions, as classical solutions. It is interesting that a single equation of motion con- 
tains various physically meaningful solutions corresponding to spaces with various topologies, 
dimensions, curvatures and sizes. 

Another direction which has been pursued is to look for relations between tensor models 
and the general relativity. In the previous paper [37] , the fluctuation spectra around a specific 
type of solutions, the Gaussian solutions, to a tensor model for one- to four- dimensional tori 
are numerically studied. It was shown that the number of the low-lying low-momentum modes 
at each momentum sector agrees exactly with what is expected from the general relativity. 

It is obvious that such number counting is not enough to definitely prove the relations, 
and more detailed comparisons are required. In fact, as will be discussed in Section 13. the 
Gaussian solutions have a natural generalization to incorporate a correspondence between the 
rank-three tensor in tensor models and the metric tensor in the general relativity. Moreover, 
they are right on the trajectory of the renormalization procedure proposed in Section [STTl and 
can be expected to play significant roles in the dynamics. Therefore the Gaussian solutions 
seem to be interesting backgrounds to perform more detailed comparisons with the general 
relativity. However, the action having such Gaussian solutions considered in the previous 
paper [37] is very complicated, and it is hard to obtain fully reliable numerical results from 
the model because of the heavy requirement of machine powers. Therefore, in the following, 
a much more simpler model having the Gaussian solutions will be given. This model contains 
a fractional inverse power of a matrix, and therefore cannot be considered as a well-defined 
action for all the values of Cabc- But, as for small fluctuations around the Gaussian solutions, 
this irregularity will not make any harms as can be seen later. 

The Gaussian solutions have the following Gaussian form, 

g-^'-^ = <5^(xi-X2), (7) 

where B, (3 are positive numerical coefficients, Xi are D-dimensional continuous coordinates, 
Xi = {xl,xf, ■ ■ ■ ,xf), and (x)^ is a short-hand notation for J2^=ii^^)'^- Only the Euclidean 
signature of space is considered in this paper. The integration measure for the contraction 
of indices is defined by J d^x. As in ([7]), throughout this paper, the symbol "on a tensor is 
used to represent the solution, which is distinguished from the dynamical variable of a tensor 
model. As can be checked easily, the algebra of functions defined by ([7]), fxifx2 = Cx-i^,x2^'^ fx^, 
is commutative nonassociative, and approaches fxifx2 ^ S^{xi — X2)fxi in the limit /3 — oo. 
Noting that ([7j) has the Poincare symmetry, and that the products of functions defined by 
/a;. = 6^{x — Xi) satisfies fxifx2 = —X2)fxi on a usual continuous space, the fuzzy space 



6 



defined by ([7]) for a finite (3 can be regarded as a commutative nonassociative deformation of a 
usual Z)-dimensional flat space. This fuzzy space is essentially the same one considered in [23] 
to investigate the one-loop properties of field theory on nonassociative space. It is important 
to note that, because of the identity form of g^^'^'^^ ([7]) has a gauge fixed form which can be 
embedded into a Euclidean tensor model. 

In the actual computation of contracting indices, however, the representation of ([7]) with 
the indices of the spatial coordinates x is inconvenient, because the momentum conservation 
coming from the translational symmetry of the solution is not obvious. One may obtain an 
expression in terms of momentum indices by insertion of the identity, 

S^{x -y) = j d^pexp {ip{x - y)) , (8) 

to all the contractions of the index x. This is equivalent to the Fourier transformation Fp^ = 
^2^^D/2 exp(ipx) for a lower index of x, and its inverse transformation F^^ = ^2J)D/2 exp(— ipx) 
for an upper index. Performing this index transformation, one obtains the Gaussian solutions 
in the momentum representation as 

C'pi,P2,P3 = ^ ^^{Pi +P2+ Ps) exp [-a {pI +pl+ pi)] , 

gP^'P^ = d^'ip.+p^), (9) 

where A and a ~ are positive numerical constants. The integration measure for the contrac- 
tion is defined by j d^p. Since the transformation is not real valued, it is not to take another 
gauge in a tensor model, but is rather to transform to a technically convenient representation 
for actual computations. Therefore, on every occasion, one has to take care of how the reality 
condition is transformed. On Cp^,p2,P3) ^^is is 

where * denotes the complex conjugation. 

To construct a tensor model with the Gaussian solutions ([2]), I follow the same procedure 
as was performed in the previous paper [37], i.e. constructing first an equation of motion 
with the Gaussian solutions, and defining an action by its square. The reason behind for 
this adhoc easy way to define an action comes from the idea that the general relativity can be 
regarded as low-energy effective phenomena associated with the spontaneous breakdown of the 
local translational symmetry [38j. In their discussions, the derivation of the general relativity 
is based on the analysis of the non-linear realization of the local coordinate transformation 
rather than starting from an action. This symmetry breaking is very similar to what occurs 
in tensor models, that is, the 0{n) symmetry, the fuzzy analogue of the general coordinate 
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transformation, breaks down at classical solutions [57]. Therefore one may guess that the 
essential properties of the low-lying spectra are determined by the symmetry breakdown, 
but not by the details of an action, although the precise correspondence to the continuum 
discussions is missing presently. 

Let me compute the following two tensors. 



TT — n p^n per' P4 

P1,P2,P3 — ^pl,p4 ^p2,p5 ^P3,P6 



K,/- ^ C,,,p„p,CP''^-^'^^ = A' (^) ' exp {-3apl) - p^), (11) 

A' (^)^ exp (-^ {pI+pI+pD^ 5^{pi+P2+P3). (12) 
Since Kp^^'^ has a diagonal form, one can safely consider its fractional powej^l, 

{K'l)p/^ = A-l {^y^exp Qap?^ 5''{p^ -p,). (13) 
Therefore the tensor W^pi,p2,P3 defined by 

Wp,,P.,P3 = Cp,,P.,Ps - (i^-^)p/k^~^)p/Kif-^)p3"^^^Pi,P^,P^ (14) 
vanishes at ([9]), if 

where K, H are defined by omitting in (ITT]) and (fT2l) , 

= CacdC'^'', (16) 

Habc = Ca/Cbe^Ccf'^, (17) 

respectively. Thus I consider 

Watc = (18) 

as the equation of motion in this paper. It should be noted that the equation of motion 
contains the Gaussian solutions irrespective of the dimension D, and is therefore giving a kind 
of unified description to all the Gaussian solutions. 

The simplest action which has (fT8l) as its equation of motion can be given by 

S = WabcW'". (19) 

It should be noted that the action (fT9|) is semi-positive definite, and if a solution to (ITSl) is 
found, this is necessarily a stable solution. Therefore there will exist no negative spectra of 
quadratic fluctuations around it. 



•'■The S^{pi — P2) is the identity matrix, any fractional power of which remains the same. 
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3 Correspondence between the rank-three tensor and 
the metric 



In the previous paper [37] , it was numerically shown that the number of the low-lying spectra 
at each momentum sector around the Gaussian solutions representing one- to four- dimensional 
fuzzy tori agrees with the expectation from the general relativity. For two-dimensional tori, 
it was observed that the number of low-lying spectra is three at p = (0,0), and one at the 
other sectors. Below I will repeat the explanation of the number distribution in the general 
relativity. 

The actual degrees of freedom of metric must be evaluated modulo the general coordinate 
transformation. The infinitesimal local coordinate transformation of the metric tensoi§ 
is given by 5g^y = V ^Vy + V^v^, where v^{x) is a local translation vector. On a flat torus, the 
transformation is expressed as 

5g^lv{.p) = iPf^Vuip) + WuV^,{p) (20) 

in the momentum representation. Since the gauge transformation is null at p = (0, 0) sector, 
there remain all the components of the metric tensor as the actual degrees of freedom, the 
number of which is given by D{D + l)/2. On the other hand, ai p ^ (0, 0) sectors, this must 
be subtracted by the number of components of f^, and the number of the actual degrees of 
freedom becomes D{D — l)/2. For D = 2, these numbers are 3 and 1, respectively, and agree 
with the numerical results mentioned above. 

It is obvious that such comparison of the numbers is not enough to definitely prove the 
relation between tensor models and general relativity. Therefore the main purpose of this 
paper is to compare the details of the corresponding modes in tensor models and general 
relativity. In the followings, I will first propose a correspondence between the rank-three 
tensor of tensor models and the metric tensor in general relativity. Then, in a slow varying 
approximation, an 0(?T,)-invariant quadratic measure of tensor models will lead to the DeWitt 
supermetric [39] with its unique parameter being fixed. This supermetric will determine the 
metric fluctuations transverse to the local coordinate transformation. Then these transverse 
metric fluctuations will determine the profiles of the fluctuation modes of tensor models under 
the correspondence, which should be compared with the numerical analysis. 

^This should not be confused with g""^ in tensor models. 
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3.1 The correspondence 



A proposal of correspondence between the rank-three tensor and the metric tensor can be 
obtained by an invariant generahzation of the solution ([7]), 

Cxux2,x3 = Bg{xi)^g{x2)^g{x3)^ exp {d{xi, X2Y + d{x2, x^f + c?(x3, Xif)] , 
^-1.-2 = 5^(a;i-X2), (21) 

where g{x) = det {g^y{x)) and d{x, y) denotes the distanc^ between the two points x, where 
the infinitesimal length is defined by ds^ = g^^dx'^dx^ as usual. The integration measure in 
contraction of indices is defined by / d^x. The parameter (3 is redundant in the sense that 
it can be absorbed into g^y. But here it is left for the later convenience to represent the 
order of g^i,, namely 0{gfj,u) ~ 1 is kept. The correspondence (12T]) is only applicable in the 
vicinity of the solution ([7]), but will be enough to analyze the small fluctuations around it. The 
expression fl?Il) respects the invariance under the general coordinate transformation, that is 
an invariant tensor such as CabcC"''"^ is invariant under the general coordinate transformation. 
This is because the distance d{x, y) is invariant from its definition, and the integration over 
x, which appears in the contraction of an index, 

Ca,MC"'^''= / d^'x./^--- , (22) 



is also invariant. It is very interesting to see that, since g^^^^'^ in fl2T]) does not depend on the 
metric tensor, the general coordinate transformation is realized within Euclidean models. 



3.2 An 0(n)-invariant measure and the DeWitt supermetric 

The measure of tensor models which was used in the numerical analysis in the previous paper 
[37] and will be used throughout this paper is the quadratic 0(n)-invariant measur^, 

dsl = 6CabcSC^''. (23) 

The corresponding measure in the space of g^u{x) can be obtained by putting the correspon- 
dence (EID into (1231). 

^The distance between a pair of points may be defined by the length of the shortest path between them. 
The definition may become singular for topological reasons or large fluctuations of g^^^, but will cause no 
problems if only small fluctuations from a background are considered as in this paper. 

II The measure may generally be added by 6Cab' SCc"''^, which will shift the parameter of the DeWitt 
supermetric discussed below. 
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By using an identity 5g = gg^'^Sg^y, the shift of Cx^^x2,x3 under an infinitesimal shift of the 
metric 6g^u{x) is given by 



2,2:3 



^g^''{xi)6g^^{xi) + ^g^"" {x2)Sg^^{x2) + ^g'"'{x3)Sg^^{x3) 



(24) 



-/36 {d{xi, X2f + d{x2., X'if + (i(x3, xif) C^^ 
Putting ([2lD into ([23]), one obtains 

d^xid^X2d^X3 \/ g{xi)^J g{x2)\/ g{x^) -g^''{xi)6g^i,{xi) + -g^''{x2)Sgf,„{x2) 



1 T 2 

+ -g^''{x3)5g^uix3) - (35 {d{xi, X2Y + d{x2, xsf + d{x3, xif) 

X exp [-2(3 {d{xi,X2f + (i(x2, x^f + d{x3, xiY)) . 



(25) 



Because of the last exponential damping factor, the integration over X2, X3 is dominated by 
the region within the distance of order -^/l//? from xi. Therefore if g^v{x) varies so slowly that 
the variation can be neglected in the length scale (3, the values of g^iy{x) in the integrand 
can be well approximated by g^u{xi). This approximation will be available for the analysis 
of the low-momentum modes around a flat background, and the systematic improvement of 
this approximation could be obtained by derivative corrections with the expansion parameter 
yiT^. Thus one obtains the corresponding metric in the space of g^u{x) in this slow varying 
approximation as 



d Xid X2d X3g{xi)'2 

3 / ~ ~ ~ \ 1 2 

-g'"'{xi)Sg^^{xi) - (35 \d{xi - Xa)^ + d{x2 - x^f + d{x3 - Xif 



X exp ( —2/3 ( d{xi — X2Y + d{x2 — x^f + d{x3 — xiY 



(26) 

where d{x) is the distance defined with the metric g^y{xi), 

d{xy =g^Axi)x^x\ (27) 
After expanding the square in the second line of fl26|) and rewriting, one obtains 

d-Aow ^ j d^xid'^X2d'^X3g{xi)^ 

1 \ . . . d 



{g^'^{x^)6g,,{x,)f - ^g^'^{x,)6g,^{x,) ^ 



ogpa{xi) 



+/^^ ( ~WB ] ^9pv{xi)5gpa{xi) 



2(3 



d 



d 



dgpuixi) dgpa{xi) 
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X exp ( -2(3 ( d{x2f + d{x2 - x^f + d{x-i 

2 



16so / d^xg{x 



Q ^ 

{g>^''{x)5g,,{x)f + -g>^'^ {x)5g,,{x)5g,.{x)^^^ 



d 



,-1 



So 



j d^x^/^ [{g^''{x)5g,,{x) f + Ag^P{x)g''- {x)5g,,{x)5gp,{x)\ 



(28) 



Here, from the first to the second equations, the Gaussian integration over X2iX^ has been 
performed and a numerical factor has been incorporated into a new numerical coefficient Sq, 
and from the second to the last equations, I have used 

d 



d 



9 



99 



dg, 



-9 



pa 



Thus the final expression of dsli^^ ( 1251) has the form of the DeWitt supermetric 
parameter has been fixed. 



(29) 
and its 



3.3 Prediction of profiles of low-lying modes in D = 2 

In this subsection, I restrict the discussions to the small fiuctuations around the D = 2 fiat 
background g^j^^, = 6^1,. I will first obtain the mode of the metric tensor transverse to the gauge 
directions (1201) . and will then use the correspondence (12T1) to obtain the corresponding mode 
in the tensor model. 

Since the gauge transformation fl20|) vanishes for a vanishing momentum, all the com- 
ponents of the metric tensor survive at the zero-momentum sector. As for a non-vanishing 
momentum, it is enough to consider the momentum {p^, 0) {p^ ^ 0), because of the rotational 
symmetry of the fiat background. Then the infinitesimal gauge transformation is given by 

5gu = 2ip^vi, 

Sg22 = 0, (30) 
Sgi2 = ip^V2. 



On the other hand, from (|28|1 . the explicit form of the DeWitt supermetric in D = 2 is 
given by 

/ 5 1 \ / 5g^, \ 
dslioni = So d'^x {6gu,6g22,Sgi2) \ 1 5 6g22 . (31) 

V 8 / V ^912 J 
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Therefore the mode transverse to the gauge transformation fl5U]) is given by 

-1 



^911 ^ 




1 





5g22 


HI 


5 





^912 J 
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(32) 



Thus the transverse mode is characterized by 

6g22 5 

^912 



- 0. (33) 

The corresponding transverse mode in the tensor model can be obtained by putting (l33ll 
into the correspondence fl2T|) . One may directly compare the obtained mode with the low- 
lying modes from the numerical analysis. However, the direct comparison through Cabc is 
inconvenient, because it has rather many cornponents and its dependence on the metric tensor 
is rather complicated and is difficult to seCJ. In fact, it is much more convenient to use 
the two-tensor Kat defined in ( TTB]) . Under the correspondence (ETj) . the dependence of the 
two-tensor on the metric tensor can be obtained as 



X exp [—P [d{xi, X2Y + d(x[, X2Y + d(xi, x-^Y + d{x[, x^Y + 2d{x2, a^s)^)] 
= Kog{xi)-^g{x[)hxp{-Pd{xi,x[f) , (34) 

where I have used the slow varying approximation mentioned previously to integrate over 
X2,xs, and Kq is a numerical factor. 

For a small perturbation from the fiat background, the variation of K is given by 



( -f35g,u ( ) (xi - X2nx, - X2r 



2 "'^ V 2 

X exp (-/9(xi - X2f) , (35) 



where the values of the metric tensor are replaced with the representative value g^u{.ixi+X2) /I) 
by using the slow varying approximation, keeping the symmetry of -ft'xi,x2 under the exchange 
of the indices. Since the numerical analysis will be performed in the momentum representation, 
it is more convenient to have 5K in the momentum representation. Performing the Fourier 
transformation mentioned in Section 12.31 one obtains 

5K = F F '^'^SK 



X2 



*In fact, I have performed the dhect comparison in some cases, and have obtained the same residts. 
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Figure 1: The contour plot of 5K^™^^{q) with = 1, /? = 1. The axes are q = {q^-,q^)- 
This choice of the parameters p^,/? is just for an example, and therefore only the qualitative 
characters of the figure is relevant. 

= Kl 5g,,{p^ + p^) {Pi - V2r{Vi - P2r exp (" Y^(Pi - P2?^ , (36) 

where is a numerical factor and 5g^y{p) is the Fourier transform of 6g^^{x). 

When a fluctuation mode 6Cabc is obtained in numerical analysis, the variation of the 
two-tensor is given by 

5Kab = SCacd C^b"^ + C^acd ^Cb"'^, (37) 

where C°j^ is a classical solution. The formulas (1371) and (l36l) make comparable the fluctuation 
modes in tensor models and the metric modes in the general relativity. 

It is worthwhile to see 6K for the transverse mode (133|) at the momentum sector (p^, 0) in 
two dimensions. Putting (15^ into and noting P1+P2 = (p^, 0) must be satisfiecO for 6K 
in fl36|l to be non-vanishing, the momentum dependence of 6K is obtained as 



(38) 

where q = (g^, g^) and ko is a numerical factor. The profile of SK^^"'"'^{q) is plotted in Figure [TJ 



At the p = (0,0) sector on a fiat torus, as can be seen in f l20|) . the gauge symmetry is 
null, and all the metric components survive. The DeWitt supermetric (ISTl) implies that the 
eigenmodes will be given by the following three ones, 

Sgu = 5g22 = 0, Sgu ^ 0, (39) 

^^^^In the continuum case, this is described by a delta function, and this must be inserted into ([38]) . But this 
is abbreviated for simphcity, since the actual interest is in tori, which have discrete momenta. 
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Figure 2: The contour plots of 6K^'^^°{q) with p = (0,0), /9 = 1 for the zero momentum 
fluctuations. The axes are q = From left to right, (15^ . (HU]1 and (HTl) . respectively. 

= ^^?22 ^ 0, 5gi2 = 0, (40) 
5gu = -S922 ^ 0, 5^12 = 0. (41) 

The corresponding SK^^^°{q) = SKg _g are plotted in Figure [2j 

4 Comarison with numerical analysis 

The main topic of this section is to compare the properties of the low-lying modes of the tensor 
model with the general relativity by using the correspondence ( 12T|) . Because of the simplicity 
of the fitting process, the actual comparison will be performed through 6K, fl36l) and ([3j 



Let me first review the method to obtain the spectra and the eigenmodes around a solution. 
Let me assume an action S{C) of a Euclidean tensor model has a classical solution C°j,^, 



dSiC) 



dC, 



abc 



= 0. (42) 

C=CO 



The spectra of the quadratic potential for the fluctuations around a classical solution 
depend on the normalization of the fluctuations. The natural normalization can be obtained 
from the measure of the path integral in the space of Cabc- The simplest 0{n) invariant 
measure is given by fl23|) assumed in the previous section. Since Cabc is symmetric under the 
permutations of its indices, the measure for its independent components is given by 

dsl = dCabJC^"' = Yl ^[(«' ^' c)] ^^i-M SC^''^'^''\ (43) 

{a,b,c) 

where (a, b, c) denotes a set of a, b, c independent of their order, C(^a,b,c) = Cabc, and m[(a, 6, c)] 



15 



is the multiplicity defiend by 

{1 a = b = c, 
3 a = bj^c, b = c^ a, or c = a^b, (44) 
6 otherwise. 

Therefore the normalized components can be given by 

SC{a,b,c) = A/m[(a,6, c)]5C(^a,b,c)^ (45) 
and the coefficient matrix of the quadratic potential for these components is 



j^{a,b,c),{d,e,f) 



1 d^S{C) 



2 9C(,.,)9C, 



{a,b,c)<J^(d,eJ) 



C=C° 



1 1 d^S{C) 



2 ^Jm[{a, b, c)]m[{d, e, /)] dC^a,b,c)9C{d,ej) 



(46) 



C=CO 



The spectra and the eigenmodes of the quadratic potential can be obtained by diagonalizing 
the symmetric matrix (H6l) . The spectra contain a number of zero modes, which come from 
the symmetry breaking of the original 0{n) symmetry of the Euclidean tensor model to the 
remaining symmetry of a classical solution C^. According to the idea of [38], these zero modes 
should be identified as the gauge symmetry non-linearly realized on a certain background 
C". Because the eigenmodes of a symmetric matrix like fH6l) are transverse to each other, the 
"physical" modes reside in the subspace transverse to the space of the zero modes. This justifies 
the reason why special attention is paid to the transverse mode to the general coordinate 
transformation in Section [331 

The numerical analysis is performed on a Windows XP64 workstation with two Opteron 
275 (2.2GHz, dual core each) processors and 8 GB memories. The C++ codes are compiled by 
the Intel C++ compiler 10.0 with OMP parallelization. NAG C Library Mark 8 and ACML 
3.6.0 & 4.0.0 are used for numerical routines. Mathematica 5.2 is used for analyzing the 
outputs. 



4.1 Fuzzy flat tori 

The classical solutions which will be considered in this subsection are Gaussian-like numerical 
solutions representing fuzzy flat tori analogous to the analytic Gaussian solutions ([7]) or ([9]) for 
fuzzy fiat spaces. The strategy of the numerical analysis is basically the same as the previous 
paper I assume such a solution for a fuzzy fiat torus has an 5*0(2) x 5*0(2) remaining 
symmetry in the same way as the translational symmetry of a usual fiat torus. Because of the 
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symmetry, momentum is a conserved quantity so that it is convenient to use the momentum 
representation in the numerical analysis. The momentum takes discrete values because of the 
finite size of the torus. Thus two-dimensional discrete momentum vectors are taken as the 
tensor index, 

p={p\p'), (p^ = -L,-L + l,-.-,L), (47) 

where L is a UV cutoff and is a positive integer. Because of the momentum conservation, a 
classical solution can be assumed to take a form, 

'^Lp2,P3 = ^Pi+P2+P3,mMPuP2,P3), (48) 

/-^^ = 5p,+,,,(o,o), (49) 

where ^(^1,^2,^3) are numerical coefficients defined only for Pi + P2 + Ps = (0,0), and sym- 
metric for the momentum variables. The measure for contracting indices is just the discrete 
sum Y.p- 



As was stressed in Section 12. 3[ the momentum representation is just a technically conve- 
nient representation, and to embed the solution into the Euclidean tensor model, one has to 
consider a coordinate representation. The discrete form of the Fourier transformation can be 
given by the transformation matrix, 

- 2ITT ^-^p (sifi^^) • 

1 / 2'7Ti \ 

where the coordintate x takes finite discrete values, 

x = {x\x'^), {x' = -L,-L + !,■■■ ,L) , (52) 

and Fx^ and Fp^ are inverse to each other. The relation between the tensors in the momentum 
and the coordinate representations is given by 

— p Pip P2p Parj (^1,'] 

^Xi,X2,X3 Xi J- X2 -'- X3 ^Pl,P2,P31 W'-'/ 

gX,,X2 ^ p^^X,p^^X2gPuP2^S^^^^^_ (54) 

Because of the solution can now be embedded into a Euclidean tensor model, if Cxi,x2,x3 
in fl55]) is real. This is satisfied, if 

MPuP2,P3y = M-PU -P2, -ps)- (55) 
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Figure 3: The plot of Kp^^p for a Gaussian-like solution for L = 5. The horizontal axes are 
P = iP^,P^)- 



In the actual process of searching for a solution, I have assumed two more ansatz. One is 
that A{pi,p2,P3) is real. The other is that A{pi,p2,P3) is invariant under the following three 
reflection symmetries, 



For example, A{{p\,pD, {pl,pl), {pI,pI)) = A{{~p\,pl), {-pl,pl), etc. These ansatz 

reduce the number of free variables, and simplify the process of solution search. These ansatz 
mean that a solution describes a fuzzy analogue of x S^, where the two S^'s have the same 
size and a reflection symmetry, and are transverse to each other. 

Stable classical solutions to the action S{C) in (IT^ can be found by searching for its local 
minima. To reduce the free variables, the above ansatz are put into S{C) to make S{A), and 
local minima of S{A) are numerically searched by a NAG routine. Actually various solutions 
to (fT8|) {S{A)=0) has been found. Since S{C) is a semi-positive definite function of C, these 
solutions are actually solutions to S{C) with respect to C, and all the fluctuation spectra 
around it are non- negative. 

For example, one can find solutions very similar to the Gaussian solutions ([9]). In Figure [3l 
the tensor Kp _p in flTBl) computed for a solution for L = 5 is plotted. The profile has a 
Gaussian-like form. 

Since momentum is the conserved quantity of the classical backgrounds, the basic strategy 
to obtain the fluctuation spectra and mode profiles around the solutions is to obtain the 
eigenvalues and eigenvectors of the submatrix in each momentum sector of the matrix P6|) . 
The technical details are given in Appendix |Xl 




(56) 
(57) 
(58) 
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As same as what was obtained in the previous paper ^7\, the spectra can be classified into 
three categories. One is the category of zero-modes, another is the low-lying modes, and the 
other is the other "heavy" modes. 

Let me first review the zero modes discussed in detail in the previous paper [37]. The 
zero modes come from the symmetry breaking of the 0{n) {n = (2L + 1)^) symmetry to 
50(2) X 5*0(2) by the classical solutions. Following the idea of [38j, these zero modes should 
be regarded as the gauge symmetry non- linearly realized on a background. Since the 0(n, R) 
transformation operates on Cx^,x2,x:j in the coordinate representation, the generators are char- 
acterized by 

M * = M 

-'■^-'■xi,X2 ^"xi,a;25 

M^,,^, = -M^,,^,, (59) 

where are the discrete coordinates (152|) . In the momentum representation, after the Fourier 
transformation Fp^ in ( |5T1) . one obtains the corresponding conditions as 

M * = M 

Mp,,p, = -Mp,,p,. (60) 

The number counting of these generators at each momentum sector was performed in the 
previous paper [37|. Taking into account the existence of the unbroken generators in the 
p = (0, 0) sector, one obtains the formula for the number of the zero modes at each momentum 
sector for D = 2 as 

^zerofr.) - / ^ (HLi (2L + 1 - \p^\) - even(p)) - 25,,(o,o) for \p^\ < 2L, 
^zero^pj - I Q otherwise, ^ ^ 

where even(p) denotes a function of whether p is an even vector p = {2i, 2j) with i,j integers 
or not, 

, - _ J 1 for an even vector p, , , 

[ otherwise. 

One can check that the number of very tiny eigenvalues in the order of machine errors < 10^^^ 
at each momentum sector agrees with the formula (IHT]) . 

The low-lying modes are the modes which are expected to describe the "low-energy" ef- 
fective dynamics of the system, and will be identified with the metric modes of the general 
relativity in the sequel. In Figure HJ the spectra of the quadratic fiuctuations at some low- 
momentum sectors around the L = 5 solution above are shown. As was obtained in the 
previous paper [37], there exist three low- lying modes at the p = (0,0) sector, and one at 
each non-zero momentum sector except p = (2, 1) (and (1,2), e.t.c). At this sector, some of 
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Figure 4: The spectra of the quadratic fluctuations around the L = 5 solution. The horizontal 
axis is \p\ = a/ (p^y + {p^Y- The solid line is 0.000115|p|'^. At the p = (0,0) sector, there 
exist three extra spectra other than the zero modes, although they cannot be distinguished 
from the origin in the flgure. At each \p\ ^ sector except p = (2, 1) (and (1, 2), e.t.c), there 
exists one low-lying spectrum. 

the "heavy" modes seem to overlap with the low-lying mode. This would be an accidental 
phenomenon caused by the smallness of L, and may be improved for larger L. Rather than 
going to larger L, which requires much larger machine power, a systematic renormalization- 
like procedure will be discussed in Section 15.21 and the spectral pattern of the low-lying modes 
will become more evident. The low-lying modes seem to be well on the line of a massless 
trajectory with the momentum dependence of \p\^. Any other momentum dependence such as 
|pp does not flt well with the data. 

Since the Einstein-Hilbert action J (Pxy/gR is a topological term in two dimensions, the 
lowest momentum dependence of the transverse modes will come from the curvature square 
term. Therefore the quartic momentum dependence of the trajectory is natural at least in two 
dimensions. 

The eigenvectors of each spectra determines the fluctuations SCabc of each mode. Putting 
these into fl37j) . one can compute the corresponding 6Kab- These are drawn for the three 
low-lying modes at the p = (0, 0) sector for the L = 5 solution in Figure [51 Comparing with 
Figure m the "lightest" mode can be identifled with the metric mode (15^ . which changes the 
relative angle of the two S'^'s. The second one is ( I40p . which is the change of the whole size. 
The last is (14T!) . which changes the relative sizes of the two S'^'s. 

Similarly, the SKab for the low-lying mode at the p = (1, 0) sector is shown in FigureEl One 
can see that this looks very similar to 5K^™'^'^{q) in Figured! which corresponds to the metric 
fluctuation transverse to the general coordinate transformation. From a statistical analysis, 
the parameter fltting with the formula fl36l) results in Sgii/Sg22 = —0.18 ± 0.01, Sgi2 — 0. 
This is in agreement with (133|) . 

These evidences show clear matching between the low-lying modes in the tensor model and 
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Figure 5: The contour plots of SKq^^g for the three low-lying modes at the p = (0, 0) sector 
of the L = 5 solution. The axes are q = [q^, q^). From left to right, the spectral values of the 
modes are 10~^^, which is the order of the machine errors, 4 x 10~^° and 8 x 10~^°, respectively. 
The figures are in good qualitative agreement with Figure [21 which correspond to the metric 
fluctuations (l39|) . (|10|) and fHTj). 



Figure 6: The contour plot of 5-^"^ _q+(i^o) for the low-lying mode at the p = (1, 0) sector of 
the L = 5 solution. The axes are q = {q^,q'^)- Its spectrum is 1.2 x 10"^. The figure agrees 
qualitatively very well with Figure [H 
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the transverse metric fluctuations in the general relativity. 

5 Renormalization procedure for tensor models 

In this section, I will discuss a Wilson's type of renormalization procedure in tensor models. 
A standard application of such a renormalization procedure is to study quantum dynamics 
of field theory by considering effective actions of renormalized fields. The study of quantum 
properties of tensor models is surely an important direction, but is out of the scope of the 
present work. In this section, I will rather use the procedure as a systematic method to single 
out slow-varying modes. 

5.1 A proposal for renormalization procedure 

The general philosophy of the Wilson's type of renormalization procedures is in the process 
of coarse graining. High frequency (momentum) modes over a cutoff are integrated out to 
produce an effective field theory of low-frequency modes below the cutoff. The non-trivial 
issue in the gravitational case is that, since the scale itself is dynamical, one cannot set a 
cutoff independently from dynamical variables. 

In a discrete system, the procedure is generally a discrete step of defining a new renor- 
malized field by averaging over some nearby site contributions. In the application to tensor 
models, since a cutoff itself must be determined from the tensor Cabc itself as a gravitational 
system, the procedure would be described by reputation of a step described by 



where i? is a function of the tensor Caba and r = 0, 1, ■ ■ ■ labels the number of the steps. 
Since the 0{n) symmetry is an essential ingredient of tensor models and it is hard to control 
its breaking, the function R must keep the invariance. One of the simplest choices of such R 
is given by the Habc in (|T6|) . 



To check that this is a kind of coarse graining process, let me consider (or ([7])) as an 
input. From (fT2l) . one can see that the Gaussian form of (Q is kept under the renormalization 
process, and obtains a discrete renormalization flow. 




(63) 




(64) 



a 



3 



(65) 
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where a^^^ denotes the parameter a in ([2]) for the r-th renormahzed tensor C*^^\ After one 
step, the value of a increases, and the fuzziness of the space becomes larger, because the 

(r) 

algebra defined by fxifx2 = CxiX2^^fxi becomes more widely spread than that determined by 
(^(^-i) Therefore the effect caused by the renormalization step can be regarded as a kind of 
coarse graining process averaging over the contributions in a nearby region. 

5.2 Classical application and trajectory of spectra 

The problem I consider in this subsection is the overlapping of the "heavy" modes on the 
low-lying modes encountered at the p = (2, 1) sector in the L = 5 spectra as in Figure |H For 
this small size of L, it is difficult to judge whether the overlapping "heavy" modes are relevant 
in the large L limit or not. In general, a renormalization procedure should improve this 
kind of ambiguity, since a renormahzed system has a larger effective size. Therefore it would 
be interesting to see how the renormalization procedure discussed in the previous subsection 
affects the spectra. 

In the path integral formulation (Euclidean), the formal strategy to obtain an effective 
action is given by 

= y"pcMe-^'^n^'^0, (66) 

where 

g-5W(c7M) ^ I p^(r-i) ^ _ ^ (c'('-i))) e-^"'-^'(^''-^0. (67) 
Therefore the effective action is obtained by 

where is the contribution from the determinant associated to the change of variables. 

In the following I will only consider the classical application of the renormalization proce- 
dure, namely, the first term of fl55]) . and leave the analysis of the full contributions for future 
study. I also consider only S^^\ which is obtained by a one step of renormalization from the 
original one. 

After one step of the classical renormalization, the equation of motion (fT8l) is given by 

W^l{C(^^)=W^,4R-'{C^'^))=0. (69) 
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Therefore a solution to wj^H = is just the renormahzation of the original solution C*^, 

(:7«o = /?(C°). (70) 
Putting fl69l) into the first derivative dW^H/dCfJr at the solution, one obtains 



ow2 



dWabc dCghi 



(71) 

C=CO 



As in Appendix \^ the spectra of the quadratic fluctuations can basically be obtained from 
the square of the matrix flTTl) with inclusion of the multiplicity factors (jH]). 

Although the above is the theoretically correct procedure for obtaining the spectra, a 
serious difficulty appears in the actual numerical computation. That is, the NAG and ACML 
routines for eigenvalues cannot produce reliable values. This can be checked by comparing the 
obtained eigenvalues with the formula fl6Tl) of the number of the zero modes. The direct reason 
for this malfunction is that the matrix dCahc/ dC^dej in f jTTI) contains very large components 
and enhances the numerical errors. The intuitive reason for such large components can be 
given as follows. The renormahzation step will pick out slow-varying important components, 
but suppress rapidly-changing unimportant ones. Therefore the derivative dCj^^/ dCdef will 
generally be very small for unimportant components of Cdef- These small values of components 
will lead to a number of large components of dCabc/dC^Jp which is the inverse of the matrix 
dcS^JdCdef. 

To overcome the problem, let me use the fact that the matrix M^^^ of the quadratic 
fluctuations can be put into a form, 

M^^) = H'^ H, (72) 

where if is a matrix and T denotes the transpose, since M^^^ is essentially a square of 0711) as 
in Appendix |X1 Even though H contains large values, the computation with H turns out to 
be much better than directly dealing with M*^^^ . To compute a few of the smallest eigenvalues 
over the zero modes, successive minimum searches of 

(££xf^ (73) 

X j_ x± 

have been performed. Here x± is the variable real vector, and varies in the vector space 
transverse to the zero modes discussed in Section 14. H when the lowest non-zero eigenvalue is 
searched. When the second lowest is searched, the x± varies in the vector space transverse to 
the lowest non-zero mode as well. This iterative procedure can be continued until a satisfactory 
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Figure 7: In the left, the two lowest non-zero spectra at each momentum sector for the 
L = 5 solution after one step of the renormalization are shown. The horizontal axis is \p\ = 
a/ (p^)'^ + {p^Y- In the right, the lowest spectra are fitted with 9.4 x 10~^|p|^. The departure 
of the spectra from the fitting line at \p\ > 4 would be reasonable on account of the existence 
of the cutoff at L = 5. 

number of the lowest eigenvalues are obtained. The inner product comes from the measure 

In Figure [3, the lowest four and two non-zero spectra at the vanishing and non- vanishing 
momentum sectors, resepctively, are plotted. The lowest trajectory can be well fitted with 
as previously. It is interesting to see that there seems to also exist a massive trajectory. The 
physical interpretation of the second trajectory is yet to be done. 



6 Summary, conclusions and discussions 



In the previous paper [37], the number distribution of the low-lying states around the Gaus- 
sian solutions representing one- to four-dimensional fuzzy fiat tori in a tensor model was 
studied, and was shown to agree with the number distribution of the metric fluctuation modes 
transverse to the general coordinate transformation in the general relativity. Although the 
agreement of the numbers for various dimensional tori is non-trivial, it is certainly not enough 
for a definite conclusion. A disadvantage in the previous numerical study was the complex- 
ity of the action, which made it difficult to go beyond the number counting. Therefore, in 
this paper, I have used a much simpler action having the Gaussian solutions, and have per- 
formed a detailed study of the profiles of the low-lying modes around the Gaussian solutions 
representing two-dimensional fuzzy flat tori. I have obtained very good matching with the 
metric fluctuation modes. The key issues in the comparison have been the proposal of the 
correspondence between the rank-three tensor in tensor models and the metric tensor in the 
general relativity, and the DeWitt supermetric derived from the invariant metric in tensor 
models under the correspondence. It has also been observed that the low-lying modes are on 
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a massless trajectory with the dependence of which would be consistent with the fact 
in the general relativity that the lowest momentum dependence will come from the R^-term 
in two dimensions. The existence of such a well-shaped trajectory also implies the existence 
of a low-momentum effective smooth manifold. These results seem to lead to the definite 
conclusion that the low-momentum effective dynamics of the small fluctuations around the 
Gaussian solutions of tensor models is described by the general relativity. 

This conclusion insists that the symmetry 0{n) is enough as the symmetry of tensor models 
to obtain the general relativity. This is rather unexpected, since the GL{n, R) symmetry is 
more directly related to the general coordinate transformation [31] . Probably, restricting only 
to the low-momentum dynamics, the distinction between the two symmetries may not be 
relevant. As future study, this can be checked by investigating the properties of tensor models 
with dynamical g"'''. 

So far the study has been concentrated on the classical aspects, but more new interesting 
results are expected to come from the study of the quantum and thermodynamic properties, as 
in the Monte Carlo analysis of noncommutative field theories |10]. A promising technical tool 
would be the renormalization procedure introduced in Section [51 Even in the one-step classical 
application, the procedure has made the low-lying spectral patterns extensively simplified and 
clearer. This kind of simplifications will be expected to occur much more in quantum mechanics 
and in thermodynamics, since quantum and thermodynamical fluctuations generally average 
over detailed classical structures. Moreover, as shown in flB^ . the Gaussian solutions are right 
on the renormalization trajectories. Therefore it can be expected that the renormalization 
procedure and the Gaussian solutions play much more important roles in quantum mechanics 
and thermodynamics of tensor models than in the classical mechanics. The optimistic hope 
is that the classification of tensor models finally reduces to that of actions having Gaussian 
solutions. Then the metric tensor in the general relativity can be regarded as the collective 
dynamical variables applicable to general tensor models under the correspondence fl?T]) . 

In addition to the study of the dynamics, it would be interesting to study the notion 
of nonassociative spaces [IS]- [25] implicitly used behind, and also pursue philosophical basis 
in view of some known principles and bounds in quantum gravity [I]-[E]- It would also be 
interesting to analytically reproduce the results obtained so far by the brute-force numerical 
computations. 
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A Details of computation of spectra 

Because of the momentum conservation of the background, the computation of spectra is 
easier in the momentum representation. From fH6|) . taking into account Wabc = at a solution 
C = one obtains 



(74) 



A slightly non-trivial matter in ( 17^ is the partial derivative of the fractional power of K 
contained in W with respect to C. One has 



a 



2 , 



-^K-^V^—^{K'-^)/. (75) 

In the momentum basis, the matrix K^J' {Kj' at C = C°) is diagonal on account of the 
momentum conservation. Moreover, from the definition f|T6|) . the matrix is semi-positive 
definite and real symmetric in the coordinate representation, and actually all the diagonal 
components are positive real at the numerical solutions C = C^. Therefore one can obtain a 
fractional power of the matrix by taking the positive branch of the fractional powers of 
these positive diagonal components, X° = (ii"")™, where /,m are positive integers. Now let 
me consider its infinitesimal variation from C = C°, 

(X° + 5Xr = (^° + 5Ky. (76) 

Taking the terms in the first order of 6X and 6K of the both sides, and using the diagonal 
forms of X° and one obtains 
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Figure 8: The graphical expression of dWabc/dCdef- A three vertex denotes C, a blob denotes 
K~9^ and a box denotes the fraction in (!78|) . A line connecting these denotes a contraction of 
indices. The multiplicity factors and the index symmetrization are abbreviated for simplicity. 



where X° and are the diagonal components of and K^, respectively. If = K^, one 
can replace the fraction in fl77|) with its obvious limit, —{K^)m~^. Thus one obtains 



dC, 



cde 



dKa' 



c=c° 



dC, 



cde 



(78) 



C=CO 



The other partial derivatives with respect to C are obvious. The final expression of dW/ dC 
is graphically shown in FigureJHl where the multiplicity factors and the index symmetrization 
are abbreviated for simplicitjo. 
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